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Abstract 

Calculating the permanent of a (0, 1) matrix is a #P-complcte problem 
but there are some classes of structured matrices for which the permanent is 
calculable in polynomial time. The most well-known example is the fixed- 
jump (0, 1) circulant matrix which, using algebraic techniques, was shown 
by Mine to satisfy a constant-coefficient fixed-order recurrence relation. 

In this note we show how, by interpreting the problem as calculating the 
number of cycle-covers in a directed circulant graph, it is straightforward 
to reprove Mine's result using combinatorial methods. This is a two step 
process: the first step is to show that the cycle-covers of directed circulant 
graphs can be evaluated using a transfer matrix argument. The second is to 
show that the associated transfer matrices, while very large, actually have 
much smaller characteristic polynomials than would a-priori be expected. 

An important consequence of this new viewpoint is that, in combina- 
tion with a new recursive decomposition of circulant-graphs, it permits ex- 
tending Mine's result to calculating the permanent of the much larger class 
of circulant matrices with non-fixed (but linear) jumps. It also permits us 
to count other types of structures in circulant graphs, e.g., Hamiltonian 
Cycles. 

1 Introduction 

Definition 1 Let A = (ojj ) be an n x n matrix. Let S n be the set of permuta- 
tions of the integers [1, . . . , n]. The permanent of A is 

n 

Perm(A) = ^ Y\ a iMi) where 77 = [ 7r ( 1 )' • • • > n ( n ^- W 

If A is a (0, 1) matrix, then A can be interpreted as the adjacency matrix of 
some directed graph G and Perm(A) is the number of directed cycle-covers in 
G, where a directed cycle-cover is a collection of disjoint cycles that cover all of 
the vertices in the graph. Alternatively, A can be interpreted as the adjacency 
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matrix of a bipartite graph G, in which case Perm(A) is the number of perfect- 
matchings in G. The permanent is a classic well-studied combinatorial object 
(see the book and later survey by Minc[T4"l 117]). 

Calculating the permanent of a (0, 1) matrix is a ^P-Complete problem 
|2(J| even when A is restricted to have only 3 non-zero entries per row [B]. The 
best known algorithm for calculating a general permanent is a straightforward 
inclusion-exclusion technique due to Ryser [14] running in G(n2 n ) time and 
polynomial space. By allowing super-polynomial space, Bax and Franklin [1] 
developed a slightly faster (although still exponential) algorithm for the (0, 1) 
case. For non-exact calculation Jerrum, Sinclair and Vigoda |12j have developed 
a fully polynomial approximation scheme for approximating the permanent of 
nonnegative matrices. 

On the other hand, for certain special structured classes of matrices one 
can exactly calculate the permanent in "polynomial time" . The most studied 
example of such a class is probably the circulant matrices, which, as discussed 
in [7], can be thought of as the borderline between the easy and hard cases. 

An n x n circulant matrix A = (a^.j) (see Figures [1] (a) and (c)) is defined 
by specifying its first row; the (i + l) s ^ row is a cyclic shift i units to the right of 
the first row, i.e., aij = oi i+fe+j-j) mod n- Let P n denote the (0, 1) n X n matrix 
with Is in positions + 1), i = 1, ... ,n — 1, and (n, 1) and Os everywhere 
else. Many of the early papers on this topic express circulant matrices in the 
form 

A n = ai P* + a 2 P° 2 + ■■■ + a k P°* (2) 

where < s\ < s 2 < • • ■ < s k < n and aj = ai iSi+ i. 

The first major result on permanents of (0, 1) circulants was due to Metropo- 
lis, Stein and Stein [13]. Let k > be fixed and A n ^ = Yli=o ^ni ^ e ^ ne n x n 
circulant matrix whose first row is composed of Is in its first k columns and Os 
everywhere else. Then |13] showed that, as a function of n, Perm(^4 n! fc) satisfies 
a fixed order constant-coefficient recurrence relation in n and therefore, could 
be calculated in polynomial time in n (after a super polynomial "start-up cost" 
in k for deriving the recurrence relation). 

This result was greatly improved by Mine who showed that it was only a 
very special case of a general rule. Let < s± < s 2 < • • • < < n be any fixed 

sequence and set A n = A n (si, . . . , s k ) = P* 1 + P° 2 H h-R* fe . In [HJQI] Mine 

proved that Perm(j4 n ) always satisfies a constant-coefficient recurrence relation 
in n of order 2 Sk — 1. Mine's theorem was proven by manipulating algebraic 
properties of A n . Note, that as mentioned by Mine, this result is difficult to 
apply for large s k since, in order to derive the coefficients of the recurrence 
relation it is first necessary to evaluate Perm(j4 n ) for n < 2(2 Sk — 1) and, using 
Ryser's algorithm, this requires time. 

Later Codenotti, Resta and various coauthors improved these results in 
various ways; e.g. in [2] showing how to evaluate sparse circulant matrices of 
size < 200; in [H|5] showing that the permanents of circulants with only three Is 
per row can be evaluated in polynomial time; in [6] showing how the permanents 
of some special sparse circulants can be expressed in terms of determinants and 
are therefore solvable in polynomial time; in [2] showing that the permanents 



2 



of dense circulants are hard to calculate and in [7] that even approximating 
the permanent of an arbitrary circulant modulo a prime p is "hard" unless 
P# P = BPP. 

In this paper we return to the original problem of Mine. Our first main 
result will be to show that if circulant matrix A n (si, . . . , s&) is interpreted as the 
adjacency matrix of a directed circulant graph C n , then counting the number of 
cycle-covers of C n using a transfer matrix approach immediately reproves Mine's 
result. In addition to rederiving Mine's original result using a combinatorial 
rather than algebraic proof this new technique permits us extend the result to 
a much larger set of circulant graphs. It will also permit us to address other 
problems, e.g., counting Hamiltonian cycles in circulant graphs, which at first 
might seem unrelated. To explain, we first need to introduce some notation. 

Definition 2 See Figure Ul Let C n 1 ' S2 '"' Sk . be the n-node directed circulant 
graph with jumps S = {si, S2, ■ ■ ■ Sfc}. (Note that this definition permits negative 
Si.) Formally, 

where 

V(n) = {0,1,. ..,n-l} 

and 

E c(n) = Ui,j) : (j - i) mod n € sj. 

Note: we will assume that S contains at least one non-negative Sj since, if all the Sj 
were negative, we could multiply them by —1 and get an isomorphic graph. Also, we 
will often write C n as shorthand for C^ l 1 ' S2r "' Sk . 

Let G = (V,E) be a graph, T C E and v £ V. Define IV T (v) to be the 
indegree of v in graph (V, T) and 0T>t(v) to be the outdegree of v in (V, T). 
T C E is a cycle-cover of G if 

Vv G V, ID T (v) = 0D T (v) = 1. (3) 

Definition 3 Let S = {s\, S2, • • • Sfc} be given. Set 

CC(n) = {T C C n : T is a cycle-cover of C n } 

and 

T{n) = \CC(n)\ = No. of cycle-covers of C n . 

Note that, by the standard correspondence mentioned previously, A n (s\, . . . , Sk) 
is the adjacency matrix of Cn' S2 ' ' Sk and T(n) = Perm(A n (si, . . . , Sfc))- So, cal- 
culating T{n) is equivalent to calculating permanents of A n (si, . . . , s^). 

There is also a well-known simple correspondence between cycle covers and 
permutations. Consider the directed complete graph with all n 2 distinct edges 
on n vertices (self-loops are permitted). Now let S n be the set of n\ permutations 
on [0, . . . , n — 1]. For a fixed permutation ir S S n , the set of edges U^T (*) ^W) 
is a cycle cover. In the other direction suppose T is a cycle cover. Define 
7T by 7r(i) = j where j is the unique vertex such that £ T. Then 7r 
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Figure 1: C n '°' : Circulant matrices (a) and (c) are, respectively, the adjacency 
matrices of circulant graphs C n 1,0,2 in (b) and (d) for n = 6, 7. 

is a permutation. This is a one-one correspondence between cycle covers and 
permutations so T{n) counts the number of permutations ir G S n restricted 
such that (vr(i) — i) mod n £ S. For example, if S = {1,2,3}, the number of 

12 3 

cycle covers in the corresponding circulant graph C n ' ' is equal to the number 
of permutations ir such that yjr(i) — i) mod n^j G {1, 2, 3}. In fact, in [19l Sec 

4.7], Stanley shows that, for fixed S, the number of such permutations satisfies 
a recurrence relation, giving an alternative derivation of Mine's result for this 
special case (but without the bound on the order of the recurrence relation 
given in [HI [16] ) . 

In [9[ [10] the authors of this paper were interested in counting spanning trees 
and other structures in undirected circulant graphs. The main tool introduced 
there was a recursive decomposition of such graphs. In Section [2] we describe 
a related recursive decomposition of directed circulant graphs. Our technique 
will be to use this decomposition to show that for some constant m there is a 
m x 1 (column) vector function T(n) such that 

Vn > 2s, T(n) = 0T{n) and T(n + 1) = AT{n) (4) 

where s is a constant to be defined later (but reduces to s = for the Mine 
formulation described previously), (3 is a 1 x m constant row- vector and A is a 
constant mx m matrix. Such an A is known as a transfer-matrix see, e.g., |19| . 

Let P(x) = ^2i=oPiX l be any polynomial that annihilates A, i.e., P(A) = 0. 
Then it is easy to see that Vn > 2s, 

J>T(n + i) = [5(j2piA n+i -Af{2s) 

i=0 \i=0 / 
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= (5A n - 2s \^Tj>iA^j f(2s) 

= {3A n - 2 ~ s QT{2s) 
= 



where denotes the m x m zero matrix and a scalar; T{n) thus satisfies the 
degree-t constant coefficient recurrence relation T(n + t) = X^i=o — 
in n. By the Cayley-Hamilton theorem, the characteristic polynomial of A 
- which has degree < m - must annihilate A, so such a polynomial exists 
and T(n) satisfies a recurrence relation of at most degree m. In our notation, 
Mine's theorem is that T(n) satisfies a recurrence relation of degree 2 s — 1. 
Unfortunately, in our construction, m = 2 2s so the characteristic polynomial 
does not suffice for our purposes. Our next step will involve showing that even 
though A is of size 2 2s x 2 2s , there is a much smaller P, of degree 2 s — 1, that 
annihilates A, thus reproving Mine's theorem. We point out that this degree 
reduction of the transfer matrix (to less than the square-root of the original size) 
is, a-priori, quite unexpected, and does not occur in the undirected-circulant 
counting problems analyzed in [91 [10] . 

One interesting consequence of this new derivation is that, unlike in Mine's 
proof, to derive the recurrence relation it is no longer necessary to start by 
spending SI f^ 2 "l time calculating the first 2 s values of T(n) using Ryser's 
method. Instead one only has to calculate A, (3, the polynomial P and the first 
2 s values of T{n) which, as we will see later, can all be done in 0(s2 4s ) time, 
reducing the start-up complexity from doubly-exponential in s to singularly 
exponential. 

Another, albeit minor, consequence of this new derivation is that it can 
also handle non-(0, 1) circulants. That is, given any matrix A n of the form ([I]), 
even when the a, are not restricted to be in {0, 1} the technique shows that 
Perm(j4 n ) satisfies a recurrence relation of degree 2 s — 1. This is only a minor 
consequence, though, since working through the details of Mine's original proof 
it is possible to modify it to get the same result. 

A much more important new consequence, and a major motivation for this 
paper, is the fact that the proof can be extended to evaluate the permanents of 
non-constant (linear) jump circulant matrices, something which has not been 
addressed before. As an example Mine's technique would not permit calculating 
the permanents of A^ n (l,n,2n), something which our new method allows. To 
explain this, we generalize Definition [2] to 

Definition 4 See Figures^ (a) and (b). Letp, s, pi,P2, ■ ■ ■ ,Pk an d s i> s 2 5 • • • > s k 
be fixed integral constants with such thatVi, < pi < p. Set S = {pin + si,p2n + 
s 2) • " " iPk n + s k}- Denote the (pn + s)-node directed circulant graph with jumps 
S by 

where 

V(n) = {0,l,...,pn + s-l} 
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(a) C\^ n (n = 4) 




(0,0) (0,1) (0,3) 

(c) Cg 1 ;"' 2 " (n = 4) 




(b) Ct"- 2n (n = 5) 




(0,0) (0,1) 



(d) C^' 2 " (n = 5) 
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Figure 2: C^™' 2 ™, a non-constant jump circulant: Solid edges are L n . Dashed 
edges are Hook(n). (a) and (b) are the circulant graphs when n = 4,5. (c) 
and (d) are corresponding lattice representations of the same graphs. The bold 
solid edges on the right of (d) are New(5) = L5 — L4. The 3 vertices{4, 9, 14} 
on the right are VN(n). Note that the dashed Hook(n) edges for both n = 4, 5 
are "independent" of n. 



and 

E c(n) = : (j - i) mod (pn + s) £ Sj. 

Figure (2a) and (2b) illustrate C^"' 2 " for n = 4,5. Figure (2c) and (2d) are 
the corresponding lattice representation, which will be introduced in section 4. 

Note that A pn+S (pin + s\,p2n + S2, ■ ■ ■ ,Pkn + Sfc) is the adjacency matrix of 
C n so, counting the cycle-covers in C n is equivalent to evaluating Perm(^4 pn+S (pi?i-|- 
si,P2n+S2, ■ ■ ■ ,pfcn+Sfc)). Our method of counting the cycle covers in C n will be 
to derive a new recursive decomposition of C n (which might be of independent 
interest) and use it to show that an analogue of (j3|) holds in the non-constant 
jump case as well; thus T(n) still satisfies a constant-coefficient recurrence re- 
lation in n. For example, Table [H shows the recurrence relation for the number 
of cycle covers m C 3n , C 3n , C 3n+1 and C 3n+1 

In the next section we describe the recursive decomposition of C n , for 
constant-jump circulants upon which our technique is based. In Section [3] we 
show how this permits easily reproving Mine's result for constant-jump circu- 
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^0,1,2 


T(n) = 2T(n - 1) - T(n - 3) 
initial values 9, 13, 12 for n = 4, 5, 6 


T(n) ~ n 
0= (l + V5)/2 


s-i0,n,2n—l 
°3n 

n l,n+l,2n 
°3n 


T(n) = 5T(n - 1) - 5T(n - 2) 

-5T(?i-3) + 6T(?i-4) 
initial values 17, 45, 113, 309 lor n = z, 3, 4, 5 


1 \n) ~ 3 


°3n+l 

n 2,n+l,2n+2 
°3n+l 


T(n) = AT(n - 1) + 5T(n - 2) 

-16T(n - 3) - 2T{n - 4) 
-8T(n - 5) - 6T(n - 6) 
+16T(n - 7) + 3T(n - 8) 
+4T(n- 9) +T(n- 10) 
initial values 31, 169, 523, 2401, 9351, 40401, 
167763, 714025, 3010351, 12766329 
for n = 2, 3, . . . , 11 


T(n) ~ ^(/. n 
^ = (l + V5)/2 
= 2 + ^ 



Table 1: The number of cycle-covers T(n) in directed circulant graphs with 
constant jumps Cn '°' and Cn' 1,2 , and with non-constant jumps C'3^ n,2n ~ 1 , and 
C3^ n+1,2n , and Cg^ 2 " - " 1 " 1 and C 2 ^ L 1,2n+2 as derived by the techniques of this 
paper. Note that for all pairs of graphs, the number of cycle covers for each 
of the graphs in the pair is the same. This is because the adjacency matrices 
for the two items in each pair are just linear circular shifts of each other so 
the permanents of their adjacency matrices are the same. The second item in 
each pair is in the form that we analyze. That is, for the constant case, having 
s\ = 0, and for the nonconstant case, having, Vi, Sj > s. 

lants. In Section[5]we then describe the generalization of the decomposition and 
the minor modifications to the proofs that are needed to extend our analysis to 
the non-constant circulants introduced in Definition |H Finally, in Section [5l we 
sketch generalizations and other uses of our technique; we first show how it can 
be extended to calculate permanents of non 0-1 circulants. We then describe 
how it can be used to calculate the moments of the the random variable count- 
ing the number of cycles in a random restricted permutations. We conclude by 
discussing how to extend the technique to counting the number of Hamiltonian 
cycles in directed circulants, extending the result of [21], which only worked for 
circulant graphs with two jumps. 

2 A Recursive Decomposition of Directed Circulant 
Graphs 

The main conceptual difficulty with deriving a recurrence relation for T(C n ) 
is that larger circulant graphs can not be built recursively out of smaller ones. 
The crucial observation, though, is that, there is another graph, L n , the lattice 
graph, that can be built recursively, and C n can then be constructed from L n 
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(c) /.;; : • and /. : : ' 

Figure 3: Cn ' , a constant jump circulant and its lattice equivalents. In (c) 
solid edges are L n ; dashed edges are Hook(n); bold edges in C 0,1,2 are New(n) 
for n = 7. 

through the addition of a constant number of edge 60. In [TO] the authors 
of this paper developed such a recursive decomposition for undirected circulant 
graphs as a tool for counting the number of spanning trees in such graphs. In 
what follows we develop a corresponding decomposition for directed circulants 
that will permit counting cycle-covers. 

We first show this for the restricted case in which S, the set of jumps, is 
constant (independent of n), where it is easy to visualize. In Section [5] we will 
see how to extend the decomposition to the more complicated case in which the 
set of jumps can depend linearly upon n, as described in Definition [2j 

We assume that = s\ < S2 < • • • < s& and set s = s&. Figure ([3]) shows 
two circulant graphs with constant jumps 0, 1, 2. Note that our assumption 
is without loss of generality, as we can choose any row of a circulant matrix 
to be the top one; for our assumption to be correct, we choose a row with a 
'1' in its first position. Equivalently, multiplying a circulant matrix by P n or 
P~ l doesn't change its permanent so we can normalize S\ = 0. For example, 
P~ 2 + P~ l + I, P~ l + I + P n and I + P n + P%, corresponding respectively, to 
graphs Cn 2, ~ 1,0 , Cn 1,0,1 and Cn' 1 ' 2 , all have the same permanent. 

To put this into context, this is very similar to the definition of Recursive families for 
undirected graphs [3] [18] , which were used for recursively building the Tutte polynomials of 
graphs in a class. 
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Definition 5 See Figure^ Let S = {s\, S2, ■ ■ ■ s&}, where the Sj are fixed 
integers. Define the n-node lattice grapiQ with jumps S by 

where 

E L (n) = [(i,j) :j-ies}. 

Now set 

Hook(n) = E c (n) - E L (n) 

and 

New(n) = E L (n + 1) - E L (n). 

Note that this implies 

L n+ \ = L n U New(n) and C n = L n U Hook(n). (5) 

The simple but important observation is that, when n is viewed as a label rather 
than as a number, Hook(n) and New(ro) are independent of the actual value of 
n. 

Lemma 1 

Hook(n) = [J { (n - j, s - j) : 1 < j < s} , 

ses 

New(n) = y){(n — s,n)}. 

ses 

Set s = Sfe . Now define 

L(n) = {0,...S-1}, 
R(n) = {n — s, . . . , n — 1}. 

Then 

Hook(n) C (R(n) x L{n)) (6) 
New(n) C (R{n) x {n}) U {(n, n)} (7) 

Important Note: In this section and the next we will always assume that 
n > 2s since this will guarantee that L(n) n R(n) = 0. Without this assumption 
some of our proofs would fail. Also note that the {(n, n)} term inNew(n) appears 
because G S. 



2 The reason for calling this a lattice graph will become visually obvious later in Definition 
1121 which generalizes this definition to the non-constant jump case. 
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(a) CCi (all edges) and Ti (solid edges) 




© © © ►© *© 

(b) CC2 (all edges) and T2 (solid edges) 




(0}.. © © 

(c) CC3 (all edges) and T3 (solid edges) 

12 

Figure 4: All of the figures are subsets of C 6 ' ' . Solid edges are in L n ; dashed 
edges are in Hook(n). The solid plus dashed edges comprise three different cycle 
covers CCi, * = 1, 2, 3 in C%. Removing the dashed Hook(n) edges leaves three le- 
gal covers T h i = 1,2,3, inL 6 . Note that s = 2, C(Ti) = C(T 2 ) = ((0,1), (0,1)) 
and C(T 3 ) = ((0,0), (0,0)). 

3 A New Proof of Mine's result 

Let CC be a cycle-cover of C n , i.e., Vv, IDt(v) = 0T>t(v) = 1. Then, from ([3|), 
in T = CC — Hook(n), almost all vertices v except (possibly) some of those that 
have an edge of Hook(ra) hanging off of them, have IDt(v) = 0Dt( v ) = 1- This 
motivates 

Definition 6 TC Ei^n) is a legal cover of L n if 

• \/v £ V, ID T (v) < 1 and 0D T (v) < 1. 

• Vv £ V — L(n), ID T (v) = 1. 

• G V -R(n), 0D T (v) = 1. 
Then, from ([5|) we have 
Lemma 2 

(a) //TC Ec{n) is a cycle-cover of C n , then 

T — Hook(n) is a legal-cover of L n . 

(b) //TC El(h +1) is a legal-cover of L n+ \, then 

T — New(n) is a legal-cover of L n . 

From the definition of legal covers we can classify and partition legal covers 
by the appropriate in/out degrees of their vertices in L(n), R(n). 
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(a) Ti U {(n - 2,n) 





Figure 5: n was increased from 6 to 7 and S = {(4,6)} C New(6) was added 
to the Ti of the previous figure. Note that, in L 7 , C(T\ U S) = C(T 2 U S) = 
since they are no longer legal covers. Also, C{T^ U S) = ( (0,0), (0,0) ). 

Definition 7 A is a binary r-tuple if 

A = (A(0),A(l),...,A(r-l)) where Vi,A(i) G {0,1}. 

Definition 8 (See Figure^). LetV be the set of2 2s possible binary tuple pairs 
(L,R) where each of L,R are, respectively, binary s tuples. 

LetT be alegal-cover ofL n . The classification of T will beC(T) = (L T ,R T ) G 
V where 

\/0<i< s, L T (i) = IDr(i) 

R T (i) = 0D T (n-l-i). 

IfTis not a legal-cover then we will write C(T) = 0. Finally, set 

C(n) = {T C Ei(n) : T is a legal cover of L n } 
C X (n) = {TeC(n) : C{T) = X} 
T x (n) = \C x (n)\ 

so Tx(n) is the number of legal- covers of L n with classification X. 

The main reason for introducing these definitions is that checking whether a 
legal cover T of L n can be completed to a cycle-cover of C n or to a legal cover 
in L n+ \ doesn't depend upon all of T but only upon its classification C(T). 

Lemma 3 See Figures [^] and 0. 

Let X = (L x , R x ) G V . Let T\ be a legal cover in L ni and T 2 be a legal cover 
ofL n2 , such that C(Ti) = C(T 2 ) = X. 
(a) Let S C Hook(n). Then, 

Ti U S is a cycle-cover of C ni 



11 



if and only if 

T2 U S is a cycle-cover of C n2 

(b) Let S C New(n). T/ien, 

C(Ti US) = C(T 2 US). 

77ia£ zs, either both Ti U S* and T 2 U S are not /eoa/ covers or, they are both 
legal covers and there is some X' £ V such that C(T\ D S) = C(T 2 U S) = X' 

Proof. To prove (a) recall that T U S is a cycle-cover of L n if and only if, 
\/v £ V, IDru5(t;) = QT>tus(v) = 1 or 

Vu G V, TDs(v) = 1 - IDt(w) and 0D 5 (u) = 1 - DD T (w) (8) 

From Lemma Q] and the definition of a legal cover we have that this is true if 
and only if 

Vi < s, lD s (i) = 1-L x (i), 

DD 5 (n-l-z) = l-i? x (i). 

and this is only dependent upon X and S and not upon n or any other properties 
of T. 

The proof of (b) is similar and omitted here. 

□ 

This lemma permits us, for X, X' £ V and S C Hook(n), to abuse the 
notations and write (XL)S) = X' to denote that, when C(T) = X, C(TU5) = 
X' . We will sometimes also write "X U S is a cycle cover 1 ' to denote that T U S 
is a cycle cover. 



Definition 9 For X, X 1 £ V, S C Hook(n) and S' C New(n) sei 

#c,s = 



1 if X U S is a cycle cover 
otherwise 



and 



Now set 



a X,X',S' 



1 C(X' U 5') = X 
otherwise 



and 



(3x= Yl te,s 
S'cHook(n) 

ax,x f = ^2 ax,x',s'- (9) 
S'cNew(n) 



Note that /3x and ax X' are constants that can be mechanically calculated. In 
fact ctx,x' is much simpler to calculate than it might initially appear seem since 
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Lemma 4 If ax x' S 1 = 1> i/ien |S"| = 1. 

Proof. In order for X' U S to be a legal cover S must include at least one edge 
that points to vertex n, so |5| > 1. From ([7]), all edges in New(ra) point to n. If 
IS"'! > 1, then IDx'us( n ) = > 1 and X' U S wouldn't be a legal cover. 

Thus ([9]) can be calculated by summing over |New(n)| = k values, instead of 
2 k values. 

Lemmas [2] and [3] immediately imply our first technical result. 
Lemma 5 

T(n) = PxT x (n) 
xev 

and 

T x {n + l) = ^ a X x'Tx(n). 
x'ep 

Let m = \V\ = 2 2s . Take any arbitrary ordering of V and define the 1 x m 
constant vector (3 = (f3x)xev and m x m constant matrix A = (axx')xx'eV- 
Finally, set T(n) = col(Tx(n))xev to be a m x 1 column vector. Then, Lemma 
[5] is exactly 

Vn > 2s, T(n) = /3T(n) and f(n + 1) = AT{n) 

which is equation As mentioned in the introduction, this immediately 

implies that T(n) satisfies a fixed-degree constant coefficient recurrence relation 
where the degree of the recurrence is at most the degree of any polynomial P(x) 
such that P(A) = 0. By the Cayley-Hamilton theorem, Q(A) = 0, Q{x) is the 
degree m = 2 2s characteristic polynomial Q(x) = det(/X — A). 

We will now see that it is possible to reduce this degree from 2 2s down to 
below 2 s . We will do this by showing that, given appropriate orderings of the 
classifications, A = (ax,x') will have a very special block diagonal format. In 
what follows, please refer to the worked example in Appendix [A] for illustration. 

Definition 10 A linear ordering on the classifications V will be called consis- 
tent if it is the lexicographic concatenation of linear orderings on its left and 
right components. 

More specifically, linear ordering "< " is consistent on V if there exist linear 
orderings "<l" and "<r" such that if X\ = (L±,Ri) and Xi = (£2,^2) we 
have X\ < X2 if and only if one of the following is true 

L\<lL2 or Li= L L 2 and Ri< R R 2 

Note that in the above definition it is not necessary for the ordering on the 
left component to be the same as the ordering on the right one (we will use this 
fact later in Lemma [7|). 
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Lemma 6 Let A = (ax,x')- If X £ V is ordered consistently, then 

(10) 



A 



(A ••• \ 
0^00 



: •- :_ 

\ o o ••• A j 



where A is some 2 s x 2 s matrix. That is, A = diag(A, A, . . . , A) where A 
contains 2 s copies of A on its diagonal. 

Proof. Suppose X = (L X ,R X ) and X' = (L X ',R X> ). 

Recall that a x ,x> = ZscNew(n) <*x,X',S where a x ,X',s = 1 if and only if 
C(X' US) = X, and is otherwise 0. 

Let L denote any binary s-tuple. Partition V up into 2 s sets of size 2 J , 
Vh = {X £ V : L x = L}. 

Note that, from Lemma [TJ if S C New(n), none of <S"s edges have endpoints 
in L(n). Intuitively, this is because edges in New(n) only connect vertices near 
the right side of the lattice and do not touch any vertices on the left side of the 
lattice. 

Thus, if ax,x',S = 1) then L x = L x . In particular this means that if 
a x,x',s = 1 then X, X' are both in the same partition set Vl- 

Now suppose that ax,x',s = 1- Let L be any other binary s-tuple and set 

X = (L, R x ) and X 1 = (Z, R x '). (11) 

Then, again using the fact that none of the endpoints of S are in L(n) we have 
that C(X' U S) = X if and only if C{X' U 5) = X so a x ,x> = &x,x>- 

When constructing matrix A = (ax,x')x,x'ev w e previously allowed any 
arbitrary ordering of "P. Ordering the X £ V consistently groups all of the X in a 
particular Vl consecutively. The observations above imply that A is partitioned 
into 2 s x 2 s blocks where each block is of size 2 s x 2 s . The non-diagonal blocks 
correspond to ax,x' where X, X' are in different partitions so all of the non- 
diagonal blocks are 0. On the other hand, the fact that ux,x' = a x X' f° r the 
X, X' defined in (jlip and the consistency of the ordering of the X tells us that 
all the diagonal blocks are copies of each other, i.e., we have proven (jlOl) . 

□ 

Corollary 1 There is a degree 2 s polynomial P(x) such that P(A) = 0. 

Proof. From the previous lemma, any polynomial P(x) that annihilates A also 
annihilates A. Since A is a 2 s x 2 s matrix, the Cayley-Hamilton theorem says 
that the characteristic polynomial P(x) of A, which is of degree 2 s , annihilates 
A. □ 
The original Mine result |14l 117) ) gave an order of 2 s — 1. We can derive 
this through a slightly more sophisticated decomposition of A. 
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Lemma 7 Let A = 

such that P(A) = 0. 



(ax,x')- Then there is a degree 2 s — 1 polynomial P(x) 



1, then we have just seen that L x = L x . Consider R x , 



Proof. If olx X' S 

R x ' and S. 

First recall from Lemma 2] that if ctx x' S = 1 then S contains exactly one 
edge. 

We claim that if ax x' S = 1 then X and X' must contain exactly the same 
number of '0's. Note that since L x = L x ' we only need to show that R x and 
R x have the same number of '0's. 

There are actually two cases. The first case is that S = {(n — s, n)}. In this 
case we are throwing away one vertex (n — s) which (because of the legality of 
X') had outdegree zero and adding a new vertex n which also has outdegree 0. 
So, the number of '0's in X U S is the same as the number of '0's in X' . 

The second case is that S = {(n — i,n)} where i < s.. Since X' U S is 
legal, vertex n — s must have already had outdegree one so throwing it away 
doesn't change the number of '0's. Adding the new vertex n with outdegree 
'0' increases the number of '0's by one. Adding edge (n — i, n) changes the 
outdegree of vertex n — i to one, decreasing the number of '0's by one. 

So, the number of '0's in X' U S is again the same as the number of '0's in 

X. 

Recall that we have that L x = L x . This suggests that we can re-order the 
entries of A so that all R x with the same number of 0's are grouped together 
(maintaining the fact that the ordering is consistent). Since there are s- 
tuples containing i '0's, A will become a block diagonal matrix of s + 1 blocks 
with block i having size (f). That is 



A 



( B 


• 








\ 





Bi ■ 













• 









V o 


• 





B n 


/ 



(12) 



where B; is a 



matrix. 



Let Pi(x) be the characteristic polynomial of B{. This has degree < 

Note that Bq and B n are both lxl matrices. By construction, Bq = B n = 
(1) so Po(x) = P n (x) = 1 — x, i.e., their characteristic polynomial is the same.. 

Because of the block diagonal form of B, P{x) = ni^t) 1 Bi(x) annihilates B. 
This polynomial has degree < ^I^o 1 (!) = 2 s — 1, proving the lemma. D 

Lemma [5] tells us that @ holds while Lemmas [6] and [7] tell us that matrix 
A is annihilated by polynomial P(x) of degree 2 s — 1. Combining them gives 
that T(n) satisfies a degree-(2 s — 1) constant coefficient recurrence relation. 



3.1 Deriving the Recurrence Relation 

We have just seen that T(n) satisfies a degree-(2 s — 1) constant coefficient re- 
currence relation. To actually derive the recurrence relation we must construct 
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• (i) a polynomial Q(n) that annihilates A = {ax,x'} and 

• (ii) the initial conditions T{n), n = 2s, 2s + 1, . . . , 2s + 2 s — 2. 

To construct Q(x), note from Lemma [6] that it suffices to calculate the 
characteristic polynomials Q(x) of matrix A. We must therefore first calculate 
the 2 2s ax,x' entries of A. 

Recall that ax,x' = SscNew(n) a x,x',s and, as noted in the proof of Lemma 
we know that if ax x' S' = 1> then S contains at most one edge. Since New(n) 
contains s edges, we can, with the appropriate data structures, calculate a.x,x' 
m 0{s 2 ) time. We can therefore calculate all the non-zero entries in the A in 
0(s 2 2 2s ) time. Finally, we can calculate Q(x) in (9(2 3s ) time, since it takes 
0(n 3 ) time to compute the characteristic polynomial of an n x n matrix 
Thus, we can calculate Q(x) in 0(2 3s ) time. 

To derive (ii), the initial conditions T(n), n = 2s, 2s + 1, . . . , 2s + 2 s — 2. 
suppose first that we already knew T{2s) and (3. Since T(n + 1) = AT(n) we 
can use the block structure from (jlOp to calculate T(n + 1) from Tin) in 0(2 3s ) 
time. It then takes only another 2 2s time to calculate T(n + 1) = j3T[n + 1). 
So, we can calculate all of the values T(n), n = 2 s + 1, 2 s + 1, . . . , 2s + 2 s — 2 
in 0(2 J ) time, improving upon the doubly exponential procedure implied by 
Mine's original result. 

It still remains to calculate Tx(2s) and (3x for all classifications X. 

Let X = (L,R). We want to calculate the number of legal covers in Lis 
with classification X. In a legal cover the number of 'O's in L must be equal to 
the number of 'O's in R. Let a\, ai, ■ ■ ■ , a% be the indices such that L(i) = and 
b\, &2, • • • , bi be the indices such that R(i) = 0. Define the set of i edges A = 
U}=i{(2s — 1 — bj,aj)}. Now define a new graph Gx as follows: (a) start with 
the lattice graph L,2g; (b) remove all edges entering vertices aj, j = 1, 2, . . . , i; 
(c) remove all edges leaving vertices 2s — 1 — bj, j = 1,2, ... ,i; (d) add the i 
edges in A. Then it is not difficult to see that T is a legal cover in Lis if and 
only if TU A is a cycle cover of Gx- Since every cycle cover of Gx must contain 
all edges in A there is a one-one correspondence between cycle covers in Gx 
and legal covers in Lis with classification X. We can therefore calculate Tx(2s) 
by calculating the permanent of the adjacency matrix of Gx which can be done 
in 0(s2 2s ) time using Ryser's algorithm. Calculating all entries in T(2s) then 
takes 0{s2 4s ) time. 

Finally, we must calculate all the fix- Let X = (L,R). If fix 7^ then the 
number of 'O's in L must be equal to the number of 'O's in R. As above, let 
a\, ai, ■ ■ ■ , a% be the indices such that L{%) = and b\,bi, . . . ,bi be the indices 
such that R(i) = 0. Now construct the i x i bipartite graph B as follows: 

Edge (j, k) G B if and only if (n — 1 — bj, a^) £ Hook(2s). 
It is not difficult to see that f3 x is exactly the number of complete matchings in 
B. We can therefore calculate (3x by evaluating the permanent of the adjacency 
matrix of B. This can be done in 0(s2 s ) time per entry and thus in 0(s2 3s ) 
time in total. 

Combining everything, we see that we can construct the recurrence relation 
and initial conditions using 0(s2 4s ) time. 
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4 Non-constant Jump Circulant Graphs 



We now extend the above definitions and lemmas to the case of non-constant 
circulants C n = c p ^^ Sl ' P2n+S2 '"'' Pkn+Sk as introduced in Definition [3J Note 
that if s = ap + (3 for some arbitrary integer a and integer (3 > 0, we can 



rewrite C n as 



C 



pi {n+a)+s[ ,p2(n+a)+s' 2 ,---,p k (n+a)+s' k 



where Vi, s[ = Si — ap. Thus, we may assume < s < p. 

Note that, using a similar argument to that in the previous section preceding 
Definition we may and do, without loss of generality, assume Vi, Sj > s. 

Analyzing non-constant jump circulants will require a change in the way 
that we visualize the nodes of C n ; until, now, as in Figure QJc), we visualized 
them as points on a line with the edges in Hook(?i) connecting the left and right 
endpoints of the line. In the non-constant jump case it will be convenient to 
visualize them as points on a bounded-height lattice, where Hook(n) connects 
the left and right boundaries of the lattice. We start by introducing a new 
graph: 

Definition 11 See Figures^ and\B Let p, s, p±,P2, ■ ■ ■ ,Pk and s\, s 2 , . . . , Sk be 

given non-negative integral constants such that \/i, < pi < p and Si > s > 0. 
Set S = {pin+si,p2n+S2, • • • ,pkn+Sk}- For u, v and integer n, set f(n; u, v) = 
un + v. Define 

C n = (V c (n), E c (n) 

where 



V(n) 



and 



u,v 



0<u<p-2 
0<v<n-l 



| IJ { (P - 1> «) :0<v<n + s-l} 



E c{n) = { ((u u vi), (u2,v 2 )) 



(u 1 ,v 1 ),(u 2 ,v 2 ) & Vein) and 
3i such that 

f(n;u 2 ,V2) - f(n;ui,vi) =pin + s« mod (pn + s) 



If s = we see that C n is simply a rectangular lattice with other regularly 
placed edges as in Figure [2J If s > then C n is a rectangular lattice with extra 
vertices extending out from its top row as in Figure El These extra vertices, 
which disturb the regularity of the lattice, are what will complicate our analysis. 

Directly from the definition we see C n is isomorphic to C n = Cp n r ^ Sl ' P2n+S2 ' " 
In particular, cycle-covers of C n are in 1-1 correspondence with cycle covers of 
C n so we can restrict ourselves to counting cycle covers of C n . We now introduce 
the generalization of Definition [5j 



Definition 12 Let p, s, pi,P2, ■ ■ ■ ,Pk and s±, s 2 , ■ ■ ■ , Sfc and S, f be as in Defi- 
nition[Tl[ Define the pn + s-node lattice graph with jumps S 



L, 



V(n),E L (n 
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where 



E L {n) 



Now set 



and 



((ui,vi), (u 2 ,v 2 )) 



(u 1 ,vi),(u 2 ,v 2 ) G Vc(n) and 
3i such that 

(a) f(n;u 2 ,v 2 ) - f(n;ux,Vi) 
and 

(b) u 2 — u\ = pi mod p 



Pin + Si mod (pn + s) 



Hook(ra) = E c {n) - E L (n) 
New(n) = E L (n + 1) - E L (n). 



Note that this implies 



L 



n+l 



and 



L n U New(n) 



L n U Hook(n). 



(13) 



We need the following intuitive lemma that was used implicitly in the con- 
stant jump case (but was so obvious there that it was not explicitly mentioned). 
The proof is straightforward but tedious and has therefore been moved to Ap- 
pendix [Bj 

Lemma 8 Let (u\,vi), (u 2 ,v 2 ) G V(n), and e = ((iii, v\), (u 2 , v 2 )), 

e G E L (n) 44> e G E L (n + 1). 

In the constant jump case, we were able to define L(n), R(n) such that 
all edges in Hook(n) went from R(n) to L(n). It turns out that this property 
remains in the non-constant jump case as well. However, as will be seen from 
the internals of the proof of LemmaO this property is a result of our assumption 
that Si > s for all i's. If this assumption did not hold, then some Hook(n) edges 
might go from L(n) to R(n). 

It is now straightforward to derive an analogue to Lemma Q] showing that 
Hook(n) and New(n) are independent of the actual value of n. Before doing so 
we will need one more definition: 



Definition 13 



NV{n)=V L {n+l)-V L {n). 



NV(n) will be the new vertices in Vt(n + 1). Note that we did not explicitly 
define this for fixed- jump circulant graphs since in the fixed-jump case NV(n) = 
VL{n + 1) — Vl(ti) = {n}, i.e., there was only the one new vertex at each step. 

Lemma 9 Set s = s^, and define 

L(n) = {(u,v) : < u < p- 1 and < v < s- 1} 
R(n) = {(u, v) : < u < p — 2 and n — s < v < n — 1} 
U {(p — 1, v) : n + s — s < v < n + s — 1} 
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(3,1) (3,3) (3,4) 




Figure 6: The graphs C^+i ' n • ( a ) an d (b) are two representations of n = 4; 
(c) and (d) are n = 5. (a) and (c) are drawn in traditional circulant format; 
(b) and (d) in lattice graph format. The bold edges in (d) are New(ra). In the 
lattice graph format the dashed edges are Hook(n) and the double-circled nodes 
denote L(n) on the left and R{n) on the right. In (b), L(n) and R{n) actually 
abut each other; L(n) are the double-circled nodes in the two leftmost columns; 
R(n) are the double-circled nodes in the three rightmost columns. As discussed 
in the text, all edges in Hook(n) go from R(n) to L(n). 
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Then 



Hook(n) C R(n) x L(n) 
New(n) C (R(n) x NV(n)) U (NV(n) x NV{n)) . 

The proof is straightforward but tedious and has therefore also been moved to 
Appendix [B] Figures [2] and [6] illustrate the lemma. 

In Section we described how to calculate the number of cycle-covers in 
constant-jump circulant graphs. Reviewing the proof, everything there followed 
directly as a consequence from the recursive decomposition of circulant graphs 
in ([5|) combined with the structural properties of the decomposition given in 
Lemma [TJ But, as we have just seen, non-constant jump circulants and their 
decompositions have exactly the same structural properties, given in (|13p and 
LemmaEJ Therefore, the entire proof developed in Section[3]can be rewritten to 
work for non-constant jump circulants. The equivalent definitions and lemmas 
needed in the non-constant jump case are stated below. 

Definition 14 T C Ei{n) is a legal cover of L n if 

• Vt; G V, IB T (v) < 1 and 0U T (v) < 1. 

• Vv £ V — L(n), ID T (v) = 1. 

• Vv £ V - R(n), 0D T (v) = 1. 
Lemma 10 

(a) //TC Ec{n) is a cycle-cover of C n , then 

T — Hook(n) is a legal-cover of L n . 

(b) IfT C El(u +1) is a legal-cover of L n+ \, then 

T — New(n) is a legal-cover of L n . 

The only major rewriting is required in the analogue to Definition [SJ The 
more complicated structure of the lattice graph in the non-constant jump case 
requires a more complicated function to map the indices of the L(n) and R{n) 
nodes. 

Definition 15 A is a binary r-tuple if 

A = 01(0), A(l), ...,A(r- 1)) where Vi, A(i) G {0, 1}. Let V be the set of 2 2 p~ s 

tuples (L,R) where L,R are two binary ps tuples. LetT be a legal-cover of L n . 

The classification of T will be C(T) = (L T , R T ) G V where 

V0 <i<ps, L T {i) = ID T (g L (i)) 
R T (i) = 0D T («^)) 

( [i/s\ , i mod s) 

/ ( L j n ~ 1 — (* m od s) ) < P — 1 

| {[i/s\ , n + s — 1 — (i mod s)) otherwise 



where 

9 L W = 
g R (l) = 
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Note: g L and g R are simply mappings of the indices of the L T (i) and R T (i) tuples to 
the nodes in L(n) and R(n). 

IfTis not a legal-cover then we will use the convention that C(T) = 0. Finally, 
set 

C(n) = {T C El{ti) : T is a legal cover of L n } 
Cx(n) = {Te£(n) : C(T) = X} 
T x {n) = \C x {n)\ 

so Tx{n) is the number of legal-covers of L n with classification X. 
Lemma 11 See Figures [?] to 

Let X = (L x , R x ) G V . Let T\ be a legal cover in L m and T2 be a legal cover 
ofL n2 , such that C(Ti) = C(T 2 ) = X. 
(a) Let S C Hook(n). Then, 

TiU S is a cycle- cover of C ni 
iff 

T2 U S is a cycle-cover of C n2 
(b) Let S C New(n). TTien, 

C(TiU5) = C(T 2 U5). 

T/wrf is, either both 7\ U S and T2 U 5 are not Zega/ covers or, they are both 
legal covers and there is some X' £ V such that C (Ti U S) = C (T2 U S) = X' 

Definition 16 For I,I'eP,SC Hook(n) and S' C New(n) set 

fl f 1 if X U S is a cycle cover ( 1 ifC(X'uS') 

Px > S ~ { oftewe and s ' " \ otherwise 

Now set 

Px = ^2 @ x ' s and ax ' x ' = ^ a x,x',s>- 
ScHook(n) S'cNew(n) 

Because NV(n) is no longer just the one vertex set {n}, Lemma H] has to 
be replaced by 

Lemma 12 Ifa x ,x>,s> = L #*en \S'\ = \NV(n)\ = p. 

The proof is very similar to that of Lemma [H We now continue with 

Lemma 13 

T(n) = £ /?xT x (n) 

7>(n + l) = ^ a x ,x<Tx{n). 
X'eP 
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(a) CCi (all edges) and Ti (solid edges) 




(b) CC2 (all edges) and T2 (solid edges) 




(c) CC3 (all edges) and T3 (solid edges) 



Figure 7: All of the figures are subsets of C A ^ +l ' n . Solid edges are 
in L n ; dashed edges are in Hook(n). The union of solid and dashed 
edges comprise different cycle covers CCi, i = 1,2,3 in C±£X ' 2n+ . Re- 
moving the dashed edges leaves three legal covers Tj, i = 1,2,3 in 
L n . C(Ti) = C(T 2 ) = ((1,1, 1,1, 0,0, 0,0), (0,1, 0,1, 0,1, 0,1)). C(T 3 ) = 
((0,1,1,1,0, 1,0,0), (0,0,0,1,1,1,0,1)). 
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(a) Ti U {(3, 14), (13, 20), (14, 4), (19, 9)} 




(b) T 2 U {(3, 14), (13, 20), (14, 4), (19, 9)} 




(c) T 3 U {(3, 14), (13, 20), (14, 4), (19, 9)} 



Figure 8: n was increased from 4 to 5 and S = {(3, 14), (13,20), (14,4), (19,9)} 
was added to Tj, of previous figure. Note that, C(T^ U S) = since it 
is no longer a legal cover (see vertex 13). C(T\ U 5) = C(T2 U 5) = 
((1,1,1,1,0, 0,0,0), (0,1, 0,0, 1,1,0,1)). 
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We reuse the concept of consistent ordering introduced in Definition [TUJ It 
is now straightforward to redo the steps of the proof of Lemma [6] to prove 

Lemma 14 Let A = (ax,x')- If X £ V is ordered consistently, then there 
exists an 2 ps x 2 ps matrix A such that, 

(A ••• \ 
A 

: '■■ ':_ 
\ ■■■ A ) 

i.e. A = diag(A,A, . . . , A) where A contains 2 ps copies of A on its diagonal. 

To see that this really is a tight generalization of Lemma [6] note that in the 
constant case p = 1 and s = and Lemma [T5] then says that the size of A is 2 s 
which is exactly the result in Lemma [6j 

The main difference between the constant-jump and non-constant jump case 
is that, in the constant-jump case we were able, in Lemma[7J to reduce the order 
of the recurrence relation from the size of A to one less than the size of A. This 
was done by using special structural properties of A. One of the facts that 
implicitly contributed to these properties was that the size of NV(n), i.e., the 
number of new vertices added at each step, was equal to one. This is not true 
in the non-constant jump case and we are therefore not able to extend Lemma 
[7] here. So, the best that we can get, from Lemma Q31 is that the recurrence 
relation T(n) satisfies a degree-(2 ps ) polynomial, an improvement of a factor of 
2 ps over the naive solution. 

For an example of such a recurrence relation, see the second set of graphs 
in Table [H 

5 Variations and Extensions 

In this section we sketch some extensions to the result in the paper as well as 
some other uses of the transfer matrix technique presented. For clarity, the 
results are only shown for constant-jump circulants. Using the techniques of 
Section H] is straightforward to generalize the results in this section to non- 
constant jump circulants as well. 

5.1 Weighted Circulants 

Until now we have assumed that our circulants are 0-1 matrices corresponding 
to being the adjacency matrices of circulant graphs. The permanents then 
counted the number of cycle covers in the corresponding circulant graphs. An 
obvious generalization is to permit the nonzero ctj to be arbitrary values. 

In this case the matrix becomes a weighted adjacency matrix. For subsets 
T of the edges in C n let weight of T be w(T) = ri(ij)6X a M- I n ^ n ^ s case the 
permanent is the sum of the weights of all cycle covers in the corresponding 
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C n , i.e., T(n) = ^TeCC„ w (T)- We can modify our technique by changing the 
definition of T x {n) in Definition [8] to 

T x (n)= Yl w ^ 

TGC x (n) 

and the definitions of Px,s an d ctx,x',S' 111 Definition [9] to be 

_ f 10 O^ 1 ) if X U S 1 is a cycle cover 

x,,s | otherwise 

and 

f if C(X' U 5') = X 

,s I q otherwise 

With these changes the rest of the derivations and analyses remain the same 
and all of the Lemmas and proofs follow accordingly. In particular, we can show 
that the permanent still satisfies a degree(2 s — 1) recurrence relation. 

5.2 Counting Cycles in Restricted Permutations 

In Section Q] we discussed how the permanent evaluates the number of restricted 
permutations using the given jumps, i.e., T(n) also counted the number of 
permutations in 

S n (S) = {ir G S n : ir[i] — i mod (n) G S}. 

We can easily modify the transfer matrix technique to answer other questions 
about these permutations. As an example, suppose that we pick a permutation 
7T uniformly at random from S n (S) and set X = # of cycles in ir. What can be 
said about the moments of XI 

First assume that, as previously, = s\ < S2 < ■ ■ ■ < Sk- Suppose now 
that for cycle cover T G CC(n) we define #c(T) to be the number of cycles 
composing cover T and set 

TC t (n)= £ (#c(T))\ 

T£CC(n) 

That is, TCo(n) = T(n) while TC\(n) is the total number of cycles summed 
over all cycle-covers in C n . Then, again by the correspondence, we have that 
the moments of X are given by 

V*>0, E(X i ) = ^-. 

TC (n) 

The interesting point is that the transfer matrix approach introduced in this 
paper can mechanically be extended through appropriate changes to the def- 
inition of Tx(n) in Definition and the definitions of f3x,s and a Xt x' S' m 
Definition O to permit showing that for every i, TCj(n) satisfies a fixed-order 
constant coefficient recurrence relation. For given, s±, S2, ■ ■ ■ , this permits, 
for example, calculating E(X) and Var(X). 
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We should note that we are only saying that for the cost functions, TCj(n), 
the transfer matrix defined by Lemma [5] exists. Lemmas [6] and [7] will no longer 
hold, though. So the degree of the recurrence relation will be 2 2s and not 2 s — 1.. 

Another complication is that in the general case TCj(n), i > 0, we may no 
longer assume that = S\ < S2 < ■ ■ ■ < Sfc. Recall that we were allowed to make 
this assumption when calculating the permanent (i = 0) because the permanent 
was invariant under rotation of rows. This is no longer true for TCj(n), i > 0. 
As an example, consider the simple circulants C® (every vertex points to itself) 
and C\ (every vertex points to its neighbor). The adjacency matrix of the first is 
I n \ the adjacency matrix of the second P^. These T are rotationally equivalent 
to each other. In both cases there is only one cycle cover; in it is the union 
of n self loops; in C\ the directed circle. So, for C°, TC\{n) = 1 while for 
C^, TC\{n) = n and their values are different. Thus, rotationally equivalent 
circulant matrices may have different values of TC\ (n) . 

We therefore need to modify our technique to work when ^ s\ by appro- 
priately modifying the definition of classifications. The major new complication 
here is that some of the edges in Hook(n) might be going from L(n) to R(n) 
rather than from R(n) to L(n). Set 

S + = {s £ S : s > 0}, s + = max sgS + s, 

S~ = {s G S : s < 0}, s~ = max seS - \s\ (if S~ = set s~ = 0) 

Now define 

L+(n) = {0, ...s+-l}, R+(n) = {n - s + , . . . , n - 1}, 
L~(n) = {0, ...<r-l}, R{n) = {n - s~, . . . , n - 1}. 

Set s = s + + s~ and let V be the set of 2 2s tuples (L + , L_, R + , i?_) where 
L~\- , L— , R-\- , R— are, respectively, binary s + , s , s + , s tuples. 

Definition 17 T C E^n) is a legal cover of L n if 

• VtJ G V, IT> T (v) < 1 and 0D T (v) < 1. 

• Vv £ V- (L+(n) UiT(n)) , ID T («) = 1. 

• V«eF - (L~{n) U R + (n)) , 0D T (v) = 1. 

Let T be a legal-cover of L n . The classification of T will now be C{T) = 
(L%,lZ,R+,RZ) G P where 

V0 < i < s+, L%(i) = IDy(i) 

i?^(i) = 0D T (n-l-i), 

V0<i<s~, i?^(i) = ID T (n-l-i), 
LL(i) = 0D T (i). 

Not that the difference between this and the previously defined classifications 
was that previously, because sq = 0, we had L~{n) = R~{n) = 0. Given these 
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Td(n) = 3TCi(n - 1) - Td{n - 2) 
-3TCi(n-3) +TCi(n-4) 

+TCi(n-5) 
initial values 22, 42, 80, 149, 274 
for n = 4, 5, 6, 7, 8 


TCi(n) 
~ .7236n0 n 


SS ~ -7236n 


,-.0,1,2 


TCi(n) = 3Td(n - 1) - 6TCi(n - 3) 
+2Td(n - 4) + 4TCi(n - 5) 
-TCi(n-6)-rCi(n-7) 
initial values 21, 32, 56, 93, 161, 275, 475 
for n = 4, 5, . . . , 10 


TCi(n) 
~ .2764n0 n 


as ~ - 2764 - 



Table 2: TC\(n) is the number of cycles summed over all cycle covers in the 
given graph with n vertices. TCo(n) = T(n) is the number of cycle covers. In 
Table [U we saw that, in both cases, TC (n) ~ (f> n where <f> = (1 + \/E)/2. 

new definitions, we can use the same transfer matrix machinery as before to 
derive recurrence relations for the Td(n). 

As an illustration recall the results from Table [1] counting the number of 
cycle covers in Cn '°' and C n ' ' 2 . Even though these two graphs are not iso- 
morphic they had the same number of cycle-covers because the adjacency matrix 
of the second is just the adjacency matrix of the first with every row (cyclically) 
shifted over one step. Since permanents are invariant under cyclic shifts both 
matrices have the same permanent which is ~ <p n where <fi = (1 + v5)/2. 

We calculated TC\ (n) for both cases with the results given in Table [2 In 
both cases we have that TC\{n) ~ cn(j) n . This means that if a permutation 
on n items is chosen at random from the corresponding distribution then, on 
average, it will have ^S ~ cn c y c l es - It i s interesting to note that that c is 
different for the two cases. 

5.3 Hamiltonian Cycles and Other Problems 

Finally, we note that a minor modification to the transfer-matrix technique 
permits using it to show that the number of Hamiltonian Cycles in a directed 
circulant graph C n also satisfies a constant-coefficient recurrence relation in 
n. This fact was previously known for undirected circulant graphs [TO] but 
doesn't seem to have been known for directed circulants, with the exception 
of the special case of in(out)-degree 2 circulants [21] , also known as two-stripe 
circulants. 

Again, as when calculating TC\{n) in the previous subsection, we may no 
longer assume that s\ = 0. We reuse the definitions of L + (n), L~(n), R + (n), R~(n) 
introduced above and define a 

Definition 18 T C Ejj(n) is a legal tour of L n if (i) T is a Hamiltonian Cycle 
of L n or (ii) 

• V«€ 7, IDt(u) < 1 and 0D T (» < 1. 
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VvG V- (L+(n)UR-(n)), IT> T {v) = 1. 



• Vt; G V- (L-(n)U R+(n)) , OD T (v) = l. 

• T contains no cycles 

Note that if T is legal and is not Hamiltonian, then T is composed of paths in 
which (i) the start of each path is in R + (n) U L~(n) (ii) the end of each path 
is in L + (n) U R~(n) and (hi) every vertex is on exactly one path (if a vertex v 
is isolated we consider it to be lying on a zero-length path that starts and ends 
at v). The classification of T will then be the union of the (start, end) pairs 
describing the starting and ending points of each path. The number of such 
classifications is finite. Furthermore, the classification of T U S where T is a 
legal tour and S C New(n) or S C Hook(n) depends only upon the classification 
of T and the edges in S. We can therefore use the method described in this 
paper to show that the number of Hamiltonian cycles in C n satisfies a recurrence 
relation. 

We conclude by noting that there is nothing particularly special about 
Hamiltonian Cycles and that the technique will enable counting many other 
structures in directed circulant graphs as well. As an example, it is not too 
difficult to modify the method to show that the number of Eulerian Tours in 
such graphs also satisfies a constant-coefficient recurrence relation in n. 

6 Conclusion 

In this paper we showed a new derivation of Mine's result [14} I17j ) that the 
permanent of parametrized circulant matrices satisfies a recurrence relation. 
Instead of being algebraic our new technique was combinatorial. We took ad- 
vantage of the fact that permanents of 0/1 matrices count the number of di- 
rected cycle covers in the matrix associated with the graph to transform the 
problem into a counting one. We were then able to decompose circulant matri- 
ces in such a way as to allow the use of the transfer matrix method to count 
the number of cycle covers. Finally, we were able to show that the transfer 
matrix was block diagonal, with all blocks being copies of each other, reducing 
the order of the characteristic polynomial of the transfer matrix (and thus of 
the corresponding recurrence relation for the permanents). 

A benefit of this new derivation is that it easily extends to the analysis of 
non-constant (linear) jump circulants, something that the original Mine result 
could not handle. It also permits counting many other properties of circulant 
graphs, e.g., the number of Hamiltonian cycles. 
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A A Worked Example for C^ 1 ' 2 

In Sections [2] and [3] we derived that T(n), the number of cycle covers in C n ' ' , 
satisfies 

Vn > 2s, T(n) = j3f(n) and f(n + 1) = AT (n) 
where (3 = (I3x)x&v and A = (a x ,x')x,X'eV- 

For Cn' 1,2 , s = 2. Definition [8] then says that every X E V is in the form 
X = (L X ,R X ) where L X ,R X G {0, l} 2 . We can therefore represent every X 
by a four-bit binary vector in which the first two bits represent L x and the last 
two R x ; there are 16 such X G V . 

Ordering the X lexicographically we calculate that (3 is 

(1000010001100001 ), 

f (4) is 

(1000021003200001 )* 
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(where the t denotes taking the transpose), and Transfer matrix A is 



^10 
110 
10 

VJ X VJ VJ 

1 






\J \J VJ VJ 








VJ VJ VJ VJ 




\ 



VJ \J \J VJ 










10 
110 
10 

\J X VJ VJ 

1 






VJ VJ VJ VJ 








VJ VJ VJ VJ 
















10 
110 
10 
1 











^ 














10 
110 
10 
1 ) 



A 



The Lexicographic ordering is consistent so, as predicted by Lemma A is 
partitioned into 16 4 x 4 blocks where all but the diagonal blocks are and all 
of the diagonal blocks are equal to some 4x4 matrix A which in this case is 

/ 1 \ 
110 
10 
\ 1 / 

Note that the lexicographic ordering on four-bit vectors also has the property 
that 

(0,0) < (0,1) < (1,0) < (1,1). 

This means that if X\ = (Li,Xi), Xi = (Z/2,i?2) and L\ = L2 then if the 
number of '0's in Ri is less than the number of '0's in R2 then X± < Xi- This 
satisfies the conditions of the ordering used in the proof of Lemma [7J which then 
implies that A should be in the form 



A 



where Bi is a (^) x matrix. We do observe this behavior with Bq = B2 = (1) 
The characteristic polynomial of Bq and Bi is Pq(x) = 




x 



1. 



and Bx 

x — 1. The characteristic polynomial of B\ is P\(x) 
This implies that 

Q(x) = Pi(x)P (x) = {x 2 -x- l)(x - 1) 

annihilates A. 

Working through the details we can then solve to find that, for C^' 1 ' 2 , T(n) - 
2T(n - 1) - T(n - 3) with initial values T(4) = 9, T(5) = 13, and T(6) = 12. 



2x 2 + 1 
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B Proofs of Lemmas [8] and [9] 

Proof of Lemma [8} 

We only prove the part. The reverse direction can be proved by the same 
argument. 

(a) If ((ui,vi), (u 2 ,v 2 )) G Ei(n), there exist i such that 

f(n;u 2 ,v 2 )-f(n;ui,vi) = pin + Sj mod (pn + s) 

and 

u 2 — u\ = pi mod p. 

If f(n;u 2 ,v 2 ) > f(n;ui,vi), then 

/(n; u 2 , v 2 ) - f(n; ui, «i) = p^n + s« and u 2 - m = pi. 

When n is increased to n + 1, 

/(n + 1; 112,172) - /(« + l;«i,«i) 
= ii 2 (n + 1) +v 2 — ui(n + 1) — vi 
= (u 2 n + v 2 - u\n - vi) + (u 2 - tti) 
= pi(n+l) + Sj 

(b) If /(n; n 2 , v 2 ) < /(n; tti, ui), then 

pn + s + f{n; u 2 , v 2 ) - f(n; Ui, v\) = pin + s, and p + u 2 - u\ = pi. 

When n is increased to n + 1, 

p(n + 1) + s + /(n + 1; u 2 , v 2 ) - f(n + 1; ui,vi) 
= p(n + 1) + s + u 2 (n + 1) + t> 2 - + 1) - t>i 
= (pn + s + u 2 n + i> 2 — 1^1^ - in) + (f> + u 2 - ui) 
= pi(n+l) + Sj 

Therefore, in both cases, ((iti,ui), (u 2 ,v 2 )) G Ei(n + 1). Q 

Proof of Lemma [9} 

We split the proof into two parts, 
(a) Hook(n) C R(n) x L(n) : 

Let e = ((ui,vi), (u 2 ,v 2 )) be an edge in Ec{n) associated with the jump 
Pin + Si. Note that e G Hook(n) if and only if 

f{n; u 2 , v 2 ) — fin; ui,v±) = pin + Sj mod pn + s and ti 2 — u\ ^ pi mod p 

There are two cases: 

(i) /(n;«i,ui) < f{n;u 2 ,v 2 ). 

=> u 2 n + v 2 = u\n + vi+pin + Si 
=> (u 2 -u\- Pi)n = Si + vi- v 2 
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M2 - ui / pi mod p implies u 2 — u\ — pi ^ 0. On the other hand, 

— 2n < Si + v i — v 2 < 2n, 

hence u 2 — u\ — pi = ±1. 

If «2 — Ui — pi = 1, then Si A- v\ — V2 = n, so v\ > n — Si > n — s. Since, by 
definition, v\ < n — 1, we also have v 2 < Sj — 1. Furthermore, iii = U2 — Pi — 1 < 
p — 1. Hence e G -R(n) x L(n). 

If U2 — u\ — pi = —1, this implies + t>i — ^2 + n = 0. We then have 
f2 > + Sj > n + s, which is not possible since it's outside the range of V2- 

(ii) f(n; ui,vi) > f(n;u 2 ,v 2 ). 

u 2 n + V2 + pn + s = u\n + V\+ pin + Sj 
=^ (u 2 - ui + p - Pi)n = Si - s - v 2 + vi . 

Similar to the previous case: 

— 2n < Si — s — V2 + v\ < 2n, 

and U2 — u\ 7^ pi mod p implies 112 — u± + p — pi = ±1. 

If ii2 — u\ +p — pi = 1, then Si — s — V2 + V1 = n. Thus v\ > n + s — Si, and 
V2 < $i — 1. Hence e G i?(n) x L(n). 

If ii2 — ui + p — pj = —1, this implies Si — s — V2 + v\ + n = 0. We then have 
U2 = u\ + pi — p — 1 < Pi — 2 < p — 1 which implies V2 < n — 1. However this 
results in t>2 = Sj — s + n + ^i>n + Sj — s>n which is not possible since it's 
outside the range of v 2 . 

Therefore, Hook(n) C (R(n) x L(n)) . 

(b) New(n) C (R(n) x JVV(n)) U (jVV(n) x NV(n)) : 
From Lemma [8] 

New(ra) C (V(n) x JW(n)) U (NV(n) x 7(n)) U (NV(n) x W(n)). 

Let e = ((u\,vi), (u 2 ,v 2 )) G New(n) be associated with jump pj(n + 1) + Sj. 

First assume G Consider the edge e' starting with (ui,v±) 

associated with jump pin + s, in Ec(n), i.e., in circulant graph C n and not 
lattice graph L n +i. Then 

e G New(n) =^ e G" £x(ra) =^ e' G Hook(n). 

So, from part (a), (ui, v%) is in R(n). Because (u 2 , v 2 ) G NV(n), ((iti, i>i), (^2,^2)) G 
J?(n) x iVV(n). 

Now assume that (2x2, 1*2) G V(n). Consider the edge e' ending with (u 2 , v 2 ) 
associated with jump p^n + Sj in Ec{n). Again 

e G New(n) =4> e G" .Ex(n) 4e£ Hook(n). 

So, from part (a), (ui,vi) G L{n). However L(n)r\NV(n) = 0, so such a (u±,v±) 
does not exist. D 
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